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Abstract 

The aim of the present report is to study massless fermion perturbations outside five-dimensional stringy 
black holes with U(l) charges. The Dirac equation was numerically solved to obtain the time profiles for the 
evolving fermion fields, and the quasinormal frequencies at intermediate times are computed by numerical 
Prony fitting and semianalytical WKB expansion at sixth order. We also computed numerically the late- 
time power law decay factors, showing that there are in correspondence with previously reported results for 
the case of boson fields in higher dimensional odd space-times. The dependence of quasinormal frequencies 
with U(l) compactification charges are studied, and the stability of this class of higher dimensional black 
holes under fermion perturbations is established. 
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1. Introduction 

The study of test field perturbations in a black 
hole background gives the opportunity to consider 
such questions as stability of the compact object un- 
der such fluctuations, as well as the form in which 
the field relaxes at very late times. The characteris- 
tic oscillation frequencies of the black hole solution 
that carried information about the particular prop- 
erties of the space-time can also be found. This fre- 
quencies, called quasinormal frequencies, have been 
widely studied by several analytical, semi-analytical 
and numerical methods since the pioneer work of 
Regge and Wheeler 0, 0, 1, i, 0, |, ©, S 0] . 

To study test fields of fermion character is al- 
ways motivated. Fermions are universally describ- 
ing matter fields, neutrinos, gravitinos and must of 
this fields are essential for the structure of black 
hole solutions coming from string/M theories. In 
this sense, a carefull analysis of the quasinormal 
spectrum of this solutions in provides a way to fix 
the parameters of the black hole and consequently 
of string/M theories 13, U, 12, 13 1. Another im- 



ability in nature. Also fermion perturbations of the 
gravitational bulk field equations for asymptotically 
AdS string/M theory black hole solutions turn out 
to be fundamental matter fields in the boundary 
conformal field theory, providing an important class 
of objects relevant for condensed matter phy sics in a 



physics 

frame of the AdS/CFT correspondence 0, [lj, \V]\ . 

We take into account the above facts as motiva- 
tions to this paper, in which we consider the evolu- 
tion of a spin- 1/2 massless Dirac field in the space- 
time of a black hole solution in five dimensions with 
three U(l) charges [l3T |. This solution can also be 
obtained, along general lines, considering intersec- 
tion of branes and strings [10|, LljJ, ll2[ , and in this 
sense the describe a family of compact objects com- 
ing from string/M theory. 

The structure of the paper is the following. Af- 
ter the section dedicated to the introduction, we 
proceed in Section 2 to present the general line ele- 
ment considered as well as to write the fundamental 
equations to use for a study of Dirac perturbations 
in five dimensional space-times, obtaining a gen- 
eral expression for the effective potential that de- 
scribes the perturbation. Section 3 is devoted to 
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portant fact is related with the stability of this black 
hole solutions, that is in direct relation with its reli- 
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the numerical study of the evolution of the consid- 
ered perturbations, as well as to present numerical 
and semi-analytical WKB results obtained for the 
quasinormal frequencies and its dependence on the 
charges of the black hole background, including the 
limit of large angular multipole numbers. In section 
4 we present the numerical results concerning the 
relaxation of the perturbation at very late times, 
and propose an analytical form for the decay fac- 
tors. Section 5 is devoted to the conclusions. 



r M = e M a 7 a . They form a Clifford algebra in five 
dimensions, i e, they satisfy the anti-conmutation 
relations 7 a ,7 f ' = -2if b , with r/ 10 = -1. 

Under a conformal transformation of the metric 
of the form : 
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the five-spinor ?p and the Dirac operator transforms 
as [3 Eg: 



2. Fundamental equations 

In five dimensions, we can obtain black hole so- 
lutions with three independent charges by the in- 
tersection of three 2-branes at a point, or from an 
intersection of one 2-brane and a 5-brane with a 
boost, i.e, with a momentum along the common 
string [l2|, EH • Upon toroidal compactification, the 
metric reads in both cases as: 

ds 2 = A(r)dt 2 + B{r)dr 2 +C(r)dn 2 3 (1) 

where d£l 2 describes the round 3-sphere and A(r) = 
r 2jZ (l- £), B(r) = /V3 (l - C(r) = 

r 2 j-1/3 
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In the above solution if at least one of the charges 
Qi is zero, the point r = is space-time singular- 
ity covered by the event horizon of the five dimen- 
sional black hole located at tjj. Moreover the case 
in which all the charges are non-zero, correspond to 
a regular black hole with an event horizon at r = Tr 
and an inner horizon at r = 0. 

The evolution of a masless spin 1/2 fermion field 
in a five dimensional curved background is de- 
scribed by the generalization of the Dirac equation: 



= 



(3) 



where y = r M V M is the Dirac operator that acts on 
the five-spinor T M are the curved space Gamma 
matrices, and the covariant derivative is defined 
as V M = dp — |a;° b 7 Q 7b,with \x and a being tan- 
gent and space-time indixes respectively, related 
by the basis of orthonormal one forms e a = e". 
The associated conection one-forms w° = ui ab obey 
u)\ A e b — 0, and the 7° are flat space-time 
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(5) 
(6) 



For a line element in the form ds 2 = ds 2 +ds 2 , where 
ds\ = g a b(x)dx a dx b and d,S2 = 9mn(y)dy m dy n the 
Dirac operator y satisfies the direct sum decompo- 
sition 



f=f x +f, 



(7) 



For two conformally related metrics, the validity of 
massless Dirac equation in one implies the validity 
of the same equation in the other. Then, it is posi- 
ble to solve the Dirac equation in the curved space 
with the line element ([T]) by performing succesive 
conformal transformations that isolate the metric 
components that depend of the angular variables, 
and applied succesively the direct sum decompo- 
sition of Dirac operator, until obtain an equiva- 
lent problem in a spacetime of the form M 2 x ^~J 3 , 
where M 2 is two-dimensional Minkowsky spacetime 
in (£,?"*) coordinates (where r* is the tortoise coor- 
dinate defined by dr* = \J~^dr) and J~J 3 is the met- 
ric describing the 3-sphere, in which the spectrum 
of massless Dirac operator is known. This proce- 
dure is general and has been applied previously to 
four dimensional stringly black hole by one of the 
authors 14 1. 

The above method allow us to obtain, for each 
component of a Dirac spinor in the manifold M 2 
defined as: 



the following equations: 



Xt(t,r) 



AfYe = 



and 



gamma matrices related with curved-space ones by 



dt 9r* 



(8) 

(9) 
(10) 



where 



A e (r) = x /-(£+l) 



This equations can be separated to obtain: 
d 2 d d 2 Q 



dt 2 dr 2 
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The above equations gives the temporal evolu- 
tion of Dirac perturbations outside the black hole 
spacetime [2(J • As the potentials V+ and V- are su- 
persymmetric to each other in the sense considered 
by Chandrasckhar in 2l(, then Q(t,r) and Xi{ti r ) 
will develop similar time evolutions and then they 
will have the same spectra, both for scattering and 
quasi-normal. At this point it should be stressed 
that for the spinor (p\ , we have these two poten- 
tials again. In the following we will work with equa- 
tion (|12[) and we eliminate the subscript (+) for the 
effective potential, defining V{r) = V+(r) . 




Figure 1: Effective potencial for a (4+l)-stringy black hole 
with i from O(bottom) to 5(top) and Qi = Q2 = Qi = 1. 



For the stringy black hole solution considered in 
this report we have for the effective potencial V(r) 
the following expression: 



V(r) 



where Xp = £ - 
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Figure (U} shows the effective potential for dif- 
ferent multipole numbers £ in the case of stringy 
black holes with Qi = Qi = Q3 = 1. The form 
of the effective potential is similar for other values 
of compactification charges, and as we can see, this 
assures the stability of the solution under fermion 
perturbations, due to its definite positive character. 



3. Time evolution of Dirac perturbations 
and quasinormal modes 

To integrate the equation (|12[) numerically we use 
the technique developed by Gundlach, Price and 
Pullin 0. 

The obtained results in the case of massless Dirac 
fields in (4 + l)-dimensional stringy black hole back- 
ground can be observed as the time-domain pro- 
files showed in Figures (0) to ([5]). In such profiles 
r = 3rjf and the time is measured in units of black 
hole event horizon. 
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Figure 2: Logaritraic plots of the time-domain evolution of 
I = and Qi = Q2 = Q3 = 1 massless Dirac perturbations. 



As is easily seen, the time evolution of Dirac per- 
turbations in ofivc dimensional stringy black holes 
space-time follows the usual dynamics for fields in 
other black hole backgrounds. After a first transient 
stage strongly dependent on the initial conditions 
and the point where the wave profile is computed, 
we observe the characteristic exponential damping 
of the perturbations associate with the quasinormal 
ringing, followed by a so-called power law tails at 
asymptotically late times. 

To compute the quasi-normal frequencies that 
dominated at intermediate times, we assume for the 
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Figure 3: Logaritmic plots of the time-domain evolution of 
1=1 and Qi = Q2 = Q3 = 1 massless Dirac perturbations. 
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Figure 4: Logaritmic plots of the time-domain evolution of 
£ = 2 and Qi = 0, Q2 = Q3 = 1 massless Dirac perturba- 
tions. 



function Q(t,r) in equation (TT2")) the time depen- 
dence: 

Ce(t,r) = Zi{r)exp(-iutt) (16) 

Then, the function Zg{r) satisfy the Schrodinger- 
type equation: 

^ + K - V(r)] Z e (r) = (17) 

The quasinormal modes are solutions of the wave 
equation (TT2")) with the specific boundary conditions 
requiring pure out-going waves at spatial infinity 
and pure in-coming waves on the event horizon. 
Thus no waves come from infinity or the event hori- 
zon. 
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Figure 5: Logaritmic plots of the time-domain evolution of 
t = 3 and Qi = 0, Q2 = Q3 = 1 massless Dirac perturba- 
tions. 



In this report we proceed to compute the quasi- 
normal frequencies using two different methods. 
The first method uses directly the numerical data 
obtained previously, and fit this data by superposi- 
tion of damping exponents. This numerical fitting 
scheme, known as Prony method, allow us to obtain 
very accurate results for the fundamental and first 
overtones 0, H| . For higher overtones is very diffi- 
cult to be implemented, because we need to do a fit- 
ting with a great number of exponentials, and also 
we need to know very well the particular time in 
which quasinormal ringing begins, a difficult point 
to be solved in general. For this reason in the fol- 
lowing we only presents the quasinormal frequen- 
cies, using this method, for the first two overtone 
number. 

The second method that we employed is a semi- 
analytical approach to solve equation (|17[) with the 
required boundary conditions, based in a WKB- 
type approximation, that can give accurate values 
of the lowest ( that is longer lived ) quasinormal 
frequencies, and was used in several papers for the 
determination of quasinormal frequencies in a vari- 
ety of systems @, @, 0, S H El 0] ■ 

Tables 1 and 2 presents the values obtained 
for the quasinormal frequencies for some multipolc 
number I and different set of charges. We also 
present the result for a particular stringy black 
holes with Q\ = 1, Qi = Q3 = 1 in Figure ©. As 
it is observed, the sixth order WKB approach gives 
results in agreement with those obtained by fitting 
the numerical integration data using Prony tech- 
nique. As it is expected the oscillation frequency 
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Table 1: Dirac quasinormal frequencies ujtjj in (4+1)- 
stringy black hole with Ql=0.5, Q2=l, Q3=1.8 for 1 = 
to I = 4 (results from sixth order WKB approximation and 
Prony fitting of time domain data). The frequencies are 
measured in units of the black hole horizon radious ru ■ 
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Figure 6: M assies Dirac quasinormal frequencies of (4+1)- 
stringy black holes with Qi = 1,Q2 = Q3 = 1- 
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Figure 7: Dependence upon charge Qi, i = 1,2, 3 of the real 
part of the masslcss Dirac quasinormal frecuencics of stringly 
black hole with the other two charges fixed Qj = 1, j 7^ i. 
The results correspond to the first overtone for multipolar 
number from I = to I = 4. 



Table 2: Dirac quasinormal frequencies u)ru in (4+1)- 
stringy black hole with Q1=0, Q2=l, Q3=l for i = to 
i = 4 (results from sixth order WKB approximation and 
Prony fitting of time domain data). The frequencies are 
measured in units of the black hole horizon radious r^j . 

increases for higher multipole and fixed overtone 
numbers. Increasing £ the magnitude of the neg- 
ative imaginary part of the fundamental overtone 
(n = 0) increases while for higher overtones the op- 
possite situation arises. For a fixed angular number 
t, the real part of the oscillation frequencies de- 
creases as the overtone number increases, and the 
magnitude of the imaginary part increases. Then, 
modes with higher overtone numbers decay faster 
than low-lying ones. As we expected for stability, 
all quasinormal frequencies calculated in this work 
have a well defined negative imaginary part. 



Figures ([7]) and show the dependence of the 
quasinormal modes with compactification charges. 
As one of the charges increases, leaving fixed the 
other two, the real part and the absolute value of 
the imaginary part of the quasinormal frequencies 
decreases, but the rates of decreasing are different. 
Then, as the charges increases, the modes labeled 
by the same angular multipole numbers quickly be- 
comes less damped and are long lived. 

An interesting case arises when the angular mul- 
tipole number becomes very large because in this 
case the first order WKB approximation becomes 
exact and we can obtain some analytical results. 
The effective potential (TT5|) for large multipole num- 
ber can be written as 

U(r)=e 2 A(r) (18) 



-0.16 



-0.18 

3 

--0.20 



-0.22 



-0.24 




0.0 



0.5 



1.0 



1.5 



2.0 



2.5 



Qi 



Figure 8: Dependence upon charge Qi, i = 1,2,3 of the 
imaginary part of the massless Dirac quasinormal frecuencies 
of stringly black hole with the other two charges fixed Qj = 
1, j ^ i. The results correspond to the first overtone for 
multipolar number from I = to I = 4. 



where A(r) = ^Sj^ and T(r) = 1 
the first order WKB approximation gives for the 
quasinormal frequencies in this limit the result: 



Then 



(19) 



power law exponents for the massless Dirac field 
considered in this paper in the space-time corre- 
sponding to higher dimensional Schwarzschild black 
hole. As we expected, the results obtained are con- 
sistent with the power law falloff mentioned in the 
previous paragraph. 




Figure 9: Tail for I = and Qi = 0, Q 2 = Q3 = 1. The 
power-law coefficients were estimated from numerical data 
represented in the dotted line. The full red line is the possible 
analytical result. 



being r m the point in which the asymptotic effec- 
tive potential ([Tclf reach its peak. This value can 
be determined as the maximum root of the equation 



-r(r) [2f(r) + r^\ +r/(r) 
ular case of equal charges Qi 
obtain the result: 



rir(r) 
dr 



In the partic- 
= Q3 = Q we 



r H yl + Q + ^l + Q + Q 2 (20) 



4. Long Time Tails 

Another important point to study is the relax- 
ation of the perturbing fermion field outside the 
black hole 0, [24J. It is a known result that in 
higher dimensional Schwarzschild black hole neutral 
massless boson fields had a late-time behavior dom- 
inated ( for a fixed r and each multipole moment 
I ) by a factor t~^ 2t+D ~ 2 ^ for odd D-dimcnsions 
and a more rapid decay factor t-( 2e + 3D ~&) for even 
D-dimensions 



25 



To study the late-time behavior, we numerically 
fit the profile data obtained in the appropriate re- 
gion of the time domain, to extract the power law 
exponents that describe the relaxation. As a test 
of our numerical fitting scheme, we obtained the 
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Figure 10: Tail for i = 4 and Q 1 = 0, Q 2 = Q3 = 1. 
The power-law coefficients were estimated from numerical 
data represented in the dotted line. The full red line is the 
possible analytical result. 

Figures © and (flT))) show two representative ex- 
amples of the results obtained, for a five dimen- 
sional stringy black hole with one charge zero and 
the other two equal to one. In general, all our nu- 
merical results indicate that, in five dimensiona, the 
decay in the Dirac case is governed by factors of 
the form t~( 2i+3 \ This result is in complete cor- 
respondence with the general result obtained for 
bosons fields in higher dimensional space times: 



if D is the dimension of the space time then the 
late time tail is described by the power-low falloff 
^ ~ t -W+D-2) for odd D - 2 5j. it seems that the 

fermion or boson character of the test field have no 
influence in the late falloff of the perturbations also 
in higher dimensions. Then, we can conclude that, 
outside five dimensional stringy black holes, as well 
as Schwarzschild black hole, the massless Dirac field 
shows identical decay at late times. 

However, we remark at this point that this depen- 
dance is only a result consistent with our numerical 
data. A simple analytical argument to support this 
late time behaviour do not exist, in contrast to the 
case for boson fields, in which the general form of 
the effective potential is suitable to expand for large 
values of the tortoise coordinate [25[ and then ex- 
tract the above power law behavior directly from 
this asymptotic expansion. The problem related 
with the analytical determination of the decay fac- 
tors for fermion perturbations in higher dimensional 
stringy black holes remains open. 



5. Concluding remarks 

In this report we considered the evolution of 
massless Dirac test field in the space time corre- 
sponding to five dimensional black hole solutions 
coming from intersecting branes in string theories. 

After the initial transient epoch, the evolutions is 
dominated by quasinormal modes, and at late times 
by a power-low falloff. We computed the quasinor- 
mal frequencies for diferent values of compactifica- 
tion charges using two different approaches, 6th or- 
der WKB and time domain integration with Prony 
fitting of the numerical data, obtaining by both 
methods very close numerical results. The obtained 
results for the dependence of the quasinormal fre- 
quencies with charges appears to be universal for 
all values of this parameter. 

We also computed the decay factors for the late 
times relaxation of fermion perturbations in five di- 
mensional stringy black holes and show that the 
result is similar for those obtained for boson field 
in higher odd dimensional space-times. 

It should be interesting to study the case of 
more higher dimensions, when remains open if the 
fermion decay picture and very late times obey the 
same power law behaviour of boson fields. Another 
important question is the analytical study of this 
late decay factors, taking into account that for po- 
tentials typical of fermion fields do not exist a sim- 



ple analytical argument to approach this problem, 
as in the case of boson fields. 

Stringy black holes obtained by intersection of 
branes are known for dimensions up to D=9, and 
then it would be interesting to analyse the evolu- 
tion of test fields in this higher dimensional back- 
grounds. In future reports we will complete this 
studies to gain a more complete knowledge about 
the evolution of fermion as well as boson perturba- 
tions in this interesting physical systems. 
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